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Introduction: One of the highlights of the TES 
observations in the polar regions has been the identifi- 
cation of a “cryptic” region in the south where CO 2 
appears to be in the form of a solid slab rather than a 
fluffy frost [1]. While the exact mechanism(s) by 
which the cryptic region is formed are still subject of 
some debate, it appears certain that a type of rapid 
metamorphism related to the high volatility of CO 2 ice 
is involved. The high volatility of CO 2 ice under mar- 
tian conditions has several Solar System analogs (N 2 
on Triton and Pluto, SO 2 on Io), thus making the mar- 
tian cryptic region somewhat less cryptic and certainly 
non-unique among planetary objects. In an end- 
member scenario, both the formation and the spectral 
properties of the cryptic region (and of other areas in 
the seasonal caps) can be quantitatively modeled by 
considering sintering of an ensemble of quasi-spherical 
CO 2 grains [2j. This model includes the special case of 
instanteneous slab formation, which occurs when the 
grains are sufficiently small (in the submicron range) 
so that their sintering timescale is short relative to the 
deposition timescale (a situation analogous to the “sin- 
tering” of water droplets falling into a pond). 

Physics of Sintering: Originally, the idea of an- 
nealed slabs of CO 2 in the martian seasonal caps was 
proposed based on an analysis of densification time- 
scales [2]. Recently, we have also evaluated the role 
played by the non-densifying sintering mechanism 
caused by vapor transport (Kelvin effect). The main 
conclusion from this recent work is that the seasonal 
CO 2 deposits on Mars rapidly metamorphose into an 
impermeable slab regardless of the initial grain size. 
The slab forming by this mechanism is expected to 
contain quasi-spherical voids that then undergo slow 
elimination by the densifying mechanisms. This 
densification process is strongly grain-size dependent, 
which can be used to explain the persistence of both 
low- and high-emissivity areas (e.g., the cryptic region 
and Mountains of Mitchell, respectively). The pro- 
posed texture for the martian CO? deposits is consis- 
tent with both TES and other observations (in particu- 
lar, the porous texture of the slab is consistent with the 
mean density of the seasonal deposits inferred from the 
MOLA data being less than the theoretical density of 
solid C0 2 [3]) and it has important consequences for 
the modeling of the physical properties of the martian 
seasonal frost. Specifically, the radiative properties of 
the frost (e.g., albedo and emissivity) are more prop- 


erly modeled by treating radiative transfer in a slab of 
solid CO 2 containing spherical voids (and other impu- 
rities such as dust grains) rather than by the usual 
model of spherical C0 2 and dust grains in vacuo. In 
the present study, this problem is tackled by finding 
the Mie solution for a spherical particle embedded in 
an absorbing host medium [4]. The Mie solution is 
then applied in a multiple scattering code [5] to com- 
pute the radiative properties of the martian CO 2 depos- 
its. 

Application of the New RT Model: The chief ad- 
vantage of the new RT model is its connection to the 
microphysical model of the cap texture and, conse- 
quently, its predictive capability. In particular, the new 
model does not require the notion of meter-sized Mie 
boulders of solid C0 2 in order to explain the high 
emissivities in the cap spectra (e.g., in the cryptic re- 
gion) but instead relates them to low porosity. The 
strong porosity dependence of the computed emissivity 
suggests that the density evolution obtained from the 
sintering model can be coupled with the radiative 
transfer calculations to predict the evolution of emis- 
sivities. Preliminary results presented at the 6 th Mars 
Conference [6] have demonstrated the capability of the 
new model to mimic the observed evolution of the cap 
emissivity (represented as the depth of the TES 25-|im 
band, BD 25 ). 

Modeling the Shape of TES Spectra: In addition 
to generating a semi-quantitative agreement with the 
evolution of BD 25 , the new model is capable of provid- 
ing a quantitative match to the shape of the TES spec- 
tra. For this more quantitative test, it is important to 
remove the component of the spectra related to atmos- 
pheric dust, which we accomplish using an approach 
based on the emission phase function (EPF) [7]. To 
date, 105 EPF-corrected spectra of the caps have been 
generated. An example of an EPF-corrected spectrum 
and the model spectra computed using the new RT 
model are shown in Figure 1. The refractive indices for 
solid C0 2 used in the calculations are as in [8], while 
the dusty spectra have been computed using optical 
constants for palagonite [9]. Application of the new 
model to match the TES spectra can in principle lead 
to maps of the porosity, void size, and dust content for 
the polar caps. In the example shown in Figure 1, a 
fairly good match to the observed spectrum is obtained 
for a slab containing 5-|im voids at 1% porosity with a 
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fractional abundance of 1-pm dust grains of 5 x 10" 4 by 
volume (the presence of dust is responsible for the 
shift of the frequency of mimimum emissivity long- 
ward of 25 pm). Of course, the RT solution for these 
parameters is non-unique, and this will necessitate the 
development of a maximum-likehood inversion 
method utilizing a priori information. The results from 
the sintering model might in fact be used as an a priori 
constraint (i.e., in a given location, the solution for 
porosity at different times should be consistent with 
the porosity evolution predicted by the sintering 
model). In addition, several important factors ne- 
glected so far (e.g., variable cap thickness, nonuniform 
density distribution with depth, presence of a thermal 
gradient) should be included, some of which could 
further improve the agreement between the coupled 
microphysical/radiative transfer model and TES obser- 
vations. Ultimately, the new model will provide a 
powerful observationally-based tool for the modeling 
of the coupled surface-atmosphere system on Mars. 



Figure 1: Black line: Observed spectrum. Dark blue 
line: computed spectrum with 5 -pm voids, no dust. 
Light blue line: computed spectrum with 1-pm voids, 
no dust. Green line: Computed spectrum with 1-pm 
voids, 5 x 10' 4 by volume of dust. Red line: Computed 
spectrum with 5-pm voids, 5 x 10' 4 by volume of dust. 
All computed spectra assume a kinetic temperature of 
146 K, slab thickness of 1 meter, 1% porosity, and 
dust grain size of 1 pm. 
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